d 3z 2_ r 2 orbital in high-T c cuprates: Excitonic spectrum, metal-insulator phase diagram, 
optical conductivity and orbital character of doped holes 
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The single-site dynamical mean-field approximation is used to solve a model of high-Tc cuprate 
superconductors which includes both d x 2_ y 2 and d 3z 2_ r 2 orbitals on the Cu as well as the relevant 
oxygen states. Both T (with apical oxygen) and T' (without apical oxygen) crystal structures 
are considered. In both phases, inclusion of the d 3z 2_ r 2 orbital is found to broaden the range of 
stability of the charge-transfer insulating phase. For equal charge-transfer energies and interaction 
strengths, the T' phase is found to be less strongly correlated than the T phase. For both structures, 
d-d excitons are found within the charge-transfer gap. However, for all physically relevant dopings 
the Fermi surface is found to have only one sheet and the admixture of d 3z 2_ r 2 into ground state 
wave function remains negligible (< 5%). Inclusion of the extra orbitals is found not to resolve the 
discrepancy between computed and observed conductivity in the insulating state. 

PACS numbers: 74.72.-h, 71.35.-y, 71.10.Fd, 71.30.+h 



I. INTRODUCTION 

More than 25 years after their discovery^ many as- 
pects of the physics of the high-T c cuprate superconduc- 
tors remain unclear i 2 -^ For a long period, researchers at- 
tempted to discuss the physics in terms of single-band 
models, including the t-J model and the one-band Hub- 
bard model4^— While much of the low-energy physics 
can be explained by single-band models with appro- 
priately chosen parameters^— many properties of the 
cuprates and other transition metal oxides require consid- 
eration of more realistic modelsi^ii 3 - - — The importance of 
the oxygen bands was stressed early on by Emery and 
Reiterjii and their ideas were encoded in the "three- 
band" modet^ - — which retains the Cu id x 2_ y 2 and O 
%Px, v orbitals on the CuC>2 plane. Early qualitative stud- 
ies of this mode l 13 ' 23- — have been followed by recent 
quantitative studies 2 ^ - — using the dynamical mean-field 
theory (DMFT) 34 ' 35 and sometimes in conjunction with 
density functional theory calculations— - — A very recent 
paper has argued that even the low energy physics may 
reveal signatures of non-Hubbard or non-i-J physics. 39 
Although the three-band model helps us in understand- 
ing various features of cuprates, it has its limitations. For 
example, the three-band model has been shown to pro- 
vide an inadequate description of the optical absorption 
at frequencies to > 2 eVi 32 ' 33 

A natural question is whether other Cu orbitals, in par- 
ticular the Cu 3d3 2 2_ r 2, play an important role. Higher 
energy spectroscopie o 40 ' 41 have detected these states, 
which may lead in particular to excitonic states in the 
spectrum—"— An early theoretical study^£ based on the 
slave boson approximation, argued that the d 3z 2_ r 2 or- 
bitals are not just admixed into the conduction band, 
but can give rise to another sheet of the Fermi surface 
at reasonable doping levels. Variations between material 
families in the energy and mixing of the d 3z 2_ r 2 orbital 



were recently argued to affect the value of the second 
neighbor hopping, thereby explaining the material de- 
pendence of T c — The comparison of theoretical and ex- 
perimental optical absorption was argued to be at least 
partially resolved by consideration of the Cu d 3z 2_ r 2 and 
apical oxygen orbitals— 

These and many other observations motivate this pa- 
per, in which we study a six-band model which includes, 
in addition to the three bands included before, the Cu 
3d 3z 2_ r 2 orbital and (depending on crystal structure) api- 
cal oxygen 2p z orbitals above and below the Cu02 plane. 
We shall present DMFT calculation of the phase dia- 
gram, spectral functions, d-d exciton spectrum, optical 
conductivity and the effect of doping into the d 3z 2_ r 2 /p z 
complex. We also study the possible importance of api- 
cal oxygen orbitals by comparing the result of T-phase 
(with apical oxygen) and T'-phase (without apical oxy- 
gen) crystal structures— 

The remainder of the paper is organized as follows. In 
Sec. |TT] we present the model and the method we em- 
ployed. Sec. IIIII gives the numerical results of the phase 
diagram and the spectral functions. Sec. IIVI discusses 
the d-d exciton spectrum. Sec. [V] shows in-plane and o 
axis optical conductivities, and in Sec.|Vl]we discuss the 
distribution of hole doping on various orbitals and its 
consequence on Fermi surfaces. We conclude in Sec. IVIII 



II. MODEL AND METHOD 

The three-band model previously considere d 21 ' 31- — in- 
volves the Cu 3d x 2_ y 2 and O 2p xy orbitals pointing to 
the Cu site in each unit cell. In this work we consider 
the six-band model, which in addition to the Cu 3d x 2_ y 2 



and planar O 2p x ,y orbitals also includes the Cu 3c?. 



:>,z' 2 



orbital and two apical oxygen sites (above and below 
the plane) with one p z orbital on each site hybridizing 
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FIG. 1: Illustration of orbitals in the six-band model for the 
T-phase crystal structure, (a) Cu d x 2_ y 2 orbital and planar 
02 pC r orbitals which couple to it; (b) Cu d 3z 2_ r 2 orbital and 
planar and apical orbitals which couple to it. The different 
colors (red and blue online) indicate the relative sign of the 
wave function. In the T' phase the apical oxygens are absent. 



Cu 3d x 2_ y 2 orbital (which are essentially what included 
in the three-band model) and panel (b) showing orbitals 
relevant to the Cu 3d 3z 2_ r 2 orbital. 

The resulting model involves six bands, and we restrict 
attention here to paramagnetic phases, so it is not neces- 
sary to write the spin dependence of the band structure 
explicitly. The band theory part of the Hamiltonian is 
thus a 6 x 6 matrix in fc-space. We neglect the periodicity 
in the z-direction, thus the Hamiltonian only has k x and 
k y dependences. The hopping integrals between Cu and 
O are also labelled on Fig. [1] we use t pd to denote the 
hopping integral between p XtV and d x 2_ y 2, t pdz between 
p x>y and d 3z 2_ r 2 and t Pzdz between p z and d 3z 2_ r 2. Our 
previous studies of three-band models shows that the pre- 
cise value and form of the oxygen-oxygen hopping do not 
affect the results in any important way^ For definite- 
ness, here we obtain estimates for the form and mag- 
nitude of the oxygen-oxygen hopping following Ref. l2lL 
which argues that the oxygen-oxygen hopping is the re- 
sult of a virtual process involving hopping on and off the 
Cu 4s orbital. Therefore we derive the six-band model 
by applying the Lowdin downfolding procedure^ to a 
model involving the six bands considered here plus a Cu 
4s band (see Appendix for details). 

We use du to denote the d x 2_ y 2 orbital, 
d z the d 3z 2_ r 2 orbital, take the basis \ip) — 

(d\\k,dzk,Pxk,Pyk,Ptt ovc ,P h z k° w ) and write tne resulting 
band-theoretic part of the Hamiltonian as 
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where 



with the Cu 3d 3z 2_ r 2 orbital. A schematic illustration 
of orbitals retained in the six-band model is shown in 
Fig. [1] with panel (a) showing orbitals relevant to the 
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and the hybridization between Cu and O orbitals: 
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We note that a linear combination of the two apical oxy- 
gen operators decouples from the problem, however for 



the ease of calculating the c-axis conductivity (Sec. [Vj 
we leave it as it is here, explicitly keeping the two apical 
oxygen orbitals separately. 

We choose t pd = 1.6eVi2i If there is cubic symme- 
try, t pdz = l/VStpd and t Pzdz = 2/y/3t pd but in the 
T phase the Cu-0 bond length is longer along the z- 
axis than x, y-axes, resulting in a smaller value of t pdz 
and t Pzdz . We follow Ref. [5l| and use t pdz = 0.5eV, 
t Pzdz = 0.8eV, t pp = 0.6eV and t PPz = 0.4eV. These 
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values are consistent with other estimates found in the 
literature^ - — The value of t PzPz has not been consid- 
ered in Ref. [5l| but since the downfolding procedure im- 
plies that t PzPz /t PPz — t pPz /tpp (see Appendix), we set 
t PzPz = 0.27eV. We note that the effect of oxyg en-oxygen 
hopping has been studied in detail in Ref. [33| and it has 
been shown that the precise values and form of oxygen- 
oxygen hopping does not change the physics in any im- 
portant way. To model the T'-phase, in which the apical 
oxygen states are absent, we set t Pz d z = t PzP = t PzPz = 0. 

We next turn to the interaction part of the Hamil- 
tonian. When more than one Cu orbital is important, 
interactions beyond the Hubbard U must be considered. 
We adopt the standard Slater-Kanamori for m 56 i 57 for the 
interacting part of the Hamiltonian H[ nt : 

H^t = U (nd h fnd {{ ,i +»d,,t n (i,,4.) 

+ U'{n dlht n dz ^ + n dhl n dzyt ) 

+ ({/' - J){n d[l:t n dz ^ + 71^,4.71(^,4.) 

- J (d^d^d^d^ + d^d\^d u + /i.e.) 

Here we have used <fjj (<i||) as the creation (annihilation) 

operator for the planar <L c a_ y a orbital, and d\ (d z ) as the 
creation (annihilation) operator for the d% z 2_ r 2 orbital. 
All the interactions are on-site so we have not written the 
site indices explicitly. We follow the conventional choice 
of U' = U—2J which comes from symmetry arguments of 
c?-orbitals. Note that in keeping with the common prac- 
tice in modelling cuprates we do not consider interactions 
on the oxygen sites. At the parameter values we consider 
the density of holes on the oxygen sites small enough that 
these interactions are not expected to be important. 

Except in the construction of the phase diagram we will 
choose the value U — 9 e V 58 ' 59 believed to be representa- 
tive of cuprates, and set the bare p and d energies equal: 
£d = £d z and e p — e Pz (where there are p z orbitals) . We 
define the bare charge-transfer energy 

A = e p - e d . (6) 

As will be seen, a large difference in d-occupancy and 
other aspects of the physics arises from difference in in- 
plane and c-axis hopping implied by the crystal structure. 
We study J = 0, 0.5 and I eV. 

We solve the model using the single-site dynami- 
cal mean field approximatio n 34 ' 35 with the continuous- 
time quantum Monte Carlo impurity solver in its 
hybridization-expansion (CT-HYB) former— To obtain 
real-frequency information we perform analytic continua- 
tion of the imaginary-axis self-energies using the method 
of Ref. [63|. The specifics are described in Refs. I32II33I . 
The key approximation is assuming that the lattice self- 
energy is momentum-independent, E(w,fc) — » 
The self energy is a matrix in the space of orbitals. Be- 
cause the interaction is local, which involves only the d 
electrons, all entries of X except the d-d components van- 
ish. 



8 p -8 d [eV] (3-band model) 
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£p-£ d -2U+5J [eV] (6-band model) 

FIG. 2: Metal-insulator phase diagram calculated in space of 
interaction strength U and p — d energy splitting e p — Ed at 
carrier density of one hole per unit cell for the six-band model 
in the T-phase at J-values indicated (lower x-axis, solid, 
dashed and dash-dotted lines) and compared to the previously 
published 32 phase diagram for the three-band model (shaded 
area, green on-line; upper x-axis). The metal-insulator phase 
transition is first order, with a region of metastability. The 
phase diagram shows the limit of stability of the metallic 
phase, A C 2 as the left-hand lines (six-band case) or left bound- 
ary of shaded region (three-band case) and the limit of sta- 
bility of the insulating phase, A c i, as the right-hand lines 
(six-band case) or right boundary of shaded area (three-band 
case) . 



III. PHASE DIAGRAM AND SPECTRAL 
FUNCTIONS 

In this section we present the metal-insulator phase 
diagram and electron spectral functions for the six-band 
model for varying choices of Hund interaction J and com- 
pare the results to the phase diagram and spectra previ- 
ously published for the three-band model^ To facilitate 
the comparison we remove the Hartrec energy by shifting 
the s-axis of the six-band model by —2U + 5 J relative 
to the three-band model. The magnitude of the Hartree 
shift can be understood as follows. In the three-band 
model the undoped compound is the d 9 state with the 
energy e d + 2e p ; adding one electron leads to the d 10 
state with the energy 2e d + 2e p + U; the two-hole state 
nearest in energy is d 9 L with energy e d + e p . Therefore 
the physical charge-transfer energy is (note that we use 
electron notation; in some of the literature the charge- 
transfer energy is defined in hole notation, without the 
U and with e d and e p reversed) 

E (d 10 ) + E (d 9 L) - 2E (d 9 ) = U-(s p - e d ) (7) 

However, in the six-band model there is an additional 
Hartree shift arising from the id^, z 2_ r 2 orbital. In this 
case the d 9 state has energy e d +2e dz +2ep+2e Pz +3U — 5 J 
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(see, e.g. Table II of Ref. (6^); the d 10 state has energy 
2e d + 2e dz + 2e p + 2e Pz + 6U - 10 J and the two-hole 
state nearest in energy is the d 9 L state whose energy is 
e d + 2e dz + £ P + 2e Pz + 3U — 5 J. The physical charge- 
transfer energy is thus 

E (d w ) + E (d 9 L) - 2E (d 9 ) 

= 3U - 5 J - (e p - e d ) 

= U - (e p - e d - 2U + 5 J) (8) 

shifted by 2[7 — 5 J relative to the three-band model. The 
spectra presented below show that six-band models with 
the same value of e p — e d — 2U + 5 J have the same en- 
ergy splitting between the non-bonding oxygen band and 
the upper Hubbard band, and that this splitting is also 
the same as would be found in a three-band model with 
charge-transfer energy e p — e d . 

The solid, dashed and dash-dotted lines in Fig. [5] show 
the phase boundaries calculated from the T-phase six- 
band model for three different values of J. The metal- 
insulator phase transition is first order— ~— with a coex- 
istence region. A C 2, the limit of stability of the metallic 
phase, is indicated by the left-hand lines in Fig. [2] The 
limit of stability of the insulating phase is denoted by 
A c i and is indicated by the right-hand lines. Once the 
Hartree shift is removed, the Hunds coupling J is seen to 
have a minor effect on the location of the phase bound- 
ary and the width of the crossover regions, although the 
crossover region is slightly narrower for larger J. 

The shaded area (green on-line) shows previously 
published 32 results for the coexistence region of the three- 
band model: the left boundary is A C 2 and the right 
boundary is A c i. Even after the Hartree shift is removed, 
the phase boundaries are displaced significantly, and the 
coexistence regime is wider. Some of the difference in 
width arises because the three-band model could be stud- 
ied to lower temperature (0.025 eV) than the six-band 
model, but the difference is larger than the thermal ef- 
fect. While a small portion of the difference in location of 
the phase boundary arises from the difference in Hartree 
shift arising from small differences in the occupancy of 
the d x 2_ y 2 orbitals, the majority of the change is due to 
non-Hartree many-body effects. In essence, in the six- 
band model the insulating phase remains stable down to 
weaker values of the effective correlation strength than in 
the three-band model. We do not have a definitive expla- 
nation of this finding at this stage; Further clarification 
of this issue is important. 

This physics is also seen in the spectral functions, pre- 
sented in Fig. [3] for the six-band model in the undoped 
case for the T [panel (a)] and T" [panel (b)] structures 
at parameters corresponding to the charge-transfer in- 
sulator (upper panels) and paramagnetic metal (lower 
panels) phases. The results are obtained by maximum- 
entropy analytic continuation of the self energies, follow- 
ing Ref. |H. We note that analytic continuation produces 
very wide tail down from —10 eV which we do not present 
since it is subject to large uncertainties while being unim- 
portant for our discussion. 
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FIG. 3: Momentum-integrated spectral functions of the six- 
band model in the undoped case (one hole per unit cell) for 
(a) T-phase and (b) T'-phase. The Fermi energy is at zero. 
Panel (a): upper part shows the A = 21.6eV(< A C 2) result 
and the lower part shows the A = 22.8eV(= A c i) result. 
Panel (b): upper part shows the A = 21.8eV(< A C 2) result 
and the lower part shows the A = 22.8eV(= A c i) result. 
Parameters: [7 = 9 eV, J = 0, T — 0.1 eV. Panel (a) upper 
part: Ed = £d. = —26.3 eV, e p = e Pz = —4.7 eV; lower 
part Ed = £d z = —26.9 eV, e p = e Pz = —4.1 eV. Panel (b) 
upper part: Ed = Sd z = —26.2 eV, e p = —4.4 eV; lower part 
£d = s dz = -27.1 eV, e p = -4.3 eV. 

The spectra of the d x 2_ y 2 and the p x . y orbitals are 
similar to that of the three-band model with compara- 
ble parameters ^ 2 ' 33 The new features are the d 3z 2_ r 2 
and (for the T structure) the apical oxygen p z orbital. 
The spectrum of the p z orbital includes a (5-function cen- 
tered at cu = e Pz because, as noted above, one linear 
combination of the p z orbital decouples. The two side- 
bands in the p z spectrum are the bonding and antibond- 
ing portions of the orbital which couple. Inclusion of 
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oxygen-oxygen hopping between different unit cells in the 
z-direction would broaden the <5-function, however this 
effect is not important for our considerations. 

The hybridization to the p z orbitals is evident in the 
spectrum of the d^, z 2_ r 2 orbital: it has mainly a two- 
peak feature which both at around the same place as the 
side-bands of the p z orbitals, although their strengths 
arc quite different. The onset of the d 3z 2_ r 2 spectrum is 
at a lower energy (around 1 eV) than that of the d x 2_ y 2 
which is due to the fact that the lattice is distorted in the 
c-direction away from the octahedron. The lower part of 
panel (a) shows the result calculated at A = A c i. The 
ground state is marginally metallic and one can see a 
narrow quasiparticle peak appears at the Fermi energy. 

Panel (b) of Fig. [3] shows the result in the T'-phase, 
where the apical oxygen orbitals are absent. We see sim- 
ilarities in the lineshape of the d x 2_ y 2 and p x ,y orbitals. 
However the d 3z 2_ r 2 spectrum is quite different: it now 
has a single peak centered at an energy slightly above 
s p , with its onset similar to panel (a). This is a main 
change induced by absence of apical oxygen sites. The 
quasiparticle peak in the lower part of panel (b) is more 
broad than that of panel (a) indicating that in this case 
the system is less strongly correlated. Further support 
for this notion comes from the values of the imaginary 
part of Matsubara-axis self energy; here larger magni- 
tudes correspond to larger values of the effective correla- 
tion strength. We find, for example, that at a doping of 
aound x = 0.1 ImS at the lowest Matsubara frequency 
is 1.6eV for the T phase and 1.3eV for the T" phase. 

In constructing the figures we selected values of A such 
that A — A C 2 was the same for the T and T" phase calcu- 
lations. We can define the charge-transfer energy empir- 
ically as the energy difference between the non-bonding 
oxygen peak and the lowest peak in the upper Hubbard 
band and the splittings in panels (a) and (b) are seen 
to be very similar. Comparison of the upper panels of 
figures (a) and (b) shows that the gap in the T" phase 
calcuation is smaller than the gap in the T phase calcu- 
lation, indicating that for comparable paramters the T" 
phase is less strongly correlated than the T phase. Com- 
parison of the upper panels of Fig. [3fa) and (b) here to 
Fig. 2(a) of Ref. .32 shows that a separation of 7 eV leads 
to metallic behavior in the three-band model but insu- 
lating behavior in the six-band model. Examination of 
data at a distance from the phase boundary in the insu- 
lating regime shows that the insulating gap is generically 
smaller in the six-band case than it is in the three-band 
case. 



IV. d-d EXCITON SPECTRUM 

In this section we discuss the d-d exciton spectrum. 
The corresponding correlation function is defined as: 

D(t) = (t t [6(t)6H0)}) (9) 



where the operator O is either the singlet exciton opera- 
tor 

or one of the triplet exciton operators 

d\ t d zi , (11) 
i(dj t d, t -dj^), (12) 
d\^d zV (13) 

Here, the du and d z operators have the same meaning as 
in Eq. ©. 

We have used the CT-HYB procedure to measure D (t) 
[Eq. (0)] in imaginary time. We note that at J = the 
singlet and triplet do not split, as expected. However 
since we are primarily concerned with the one-hole state, 
even at non-zero J (up to 1 eV) the difference between 
the singlet and triplet is negligible. Moreover, we have 
found (not shown) that J induces very small effect on 
the exciton correlation function. We have cross-checked 
the lack of J-dependence by exactly diagonalizing an iso- 
lated CuC>6 octahedron. We therefore focus on the J = 
results only. Fig. [4] shows the results of excitonic corre- 
lation functions for the T (upper panel) and T' (lower 
panel) structures, on the imaginary-time axis on a semi- 
log scale for metallic (blue traces) and insulating (red 
traces) situations. 

D(t) is related to the real axis spectral function D(ui) 

by 

»M = / d "CT <"> 

In the regime near r = j3, D(t) drops exponentially, 
as expected if the real-axis exciton spectrum includes a 
5 function-like peak. To find the energy of the peak we 
perform a fit of D(t) to D(j3) exp(— Slj ■ r) where 5ui 
indicates the binding energy of the peak. The results are 
shown on Fig. U] with the fitting parameter 8uj indicated. 
We note that although for insulating cases the exciton 
spectra peaks inside the optical gap, in the metallic cases 
the exciton has a peak with an even larger gap, meaning 
that it moves to a slightly higher binding energy. 

The correlation function D(t) calculated from CT- 
HYB is essentially exact: it includes all quantum fluctu- 
ations. It is interesting to view this correlation function 
in diagrammatic terms as a combination of bubble dia- 
gram (convolution of interacting Green's function) and 
vertex correction. The exciton corresponds to moving an 
electron from a d$ z 2_ r 2 orbital to a d x 2_ y 2 orbital. The 
corresponding bubble diagram is 

B d (iQ n ) = -- G d|| {iu n )Gdz ( iu n + i£l n ). (15) 
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FIG. 4: Excitonic correlation functions on imaginary-time 
axis D(t) of the six-band model in the undoped case for (a) 
T-phase and (b) T'-phase. Note the semi-log scale. Two dif- 
ferent A values with a A < A C 2 (empty square with lines, red 
online) and a A = A c i (filled circle with lines, blue online) 
are shown in each panel. The parameters are the same as in 
Fig. [3] Lines without symbols: a fit to D(/3) exp(— 8uj ■ r). 
The fitted 8uj is indicated on the figure. 



FIG. 5: Exciton spectrum (dashed line, blue on-line) calcu- 
lated for six-band model at carrier concentration of one hole 
per unit cell for (a) T structure (b) and T' structure and com- 
pared to d 3z 2 _ r 2 / d x 2 _ y 2 joint density of states Bd (solid line, 
red on-line) and total joint density of states B to t (dash-dotted 
line, magenta on-line). The parameters are indicated on the 
figure and are the same as in Fig. [3] Note that the spectra 
have been rescaled to facilitate comparison of structures. 



or, on the real frequency axis, 

B d (u) = J dco'A dll (u')A^ (u + w') [/(w + J) - /(«')] 

.(16) 

where /(w) is the Fermi function. Bd(io) is the joint d- 
density of states of the d x 2_ y 2 and d 3z 2_ r 2 orbitals. 

It will be useful in our subsequent discussion to define 
the total joint density of states as 

BtoM = J dio'A tot (io')A tot (io + lo') [f(uj + u 1 ) - /(w')] 

(17) 



where A tot (w) is the total spectral function. 

In order to compare with these real frequency func- 
tions, we have analytically continued the D(t) data using 
the maximum entropy method. 64 Results are presented in 
Fig- EI The broadening of D(uj) is due to the uncertainty 
of the analytic continuation procedure but the center of 
the peak is consistent with the exponential fit shown in 
Fig. |H It is clear from Fig. [5] that in the insulating case 
the exciton spectrum has a peak inside the optical gap, 
while in the metallic case where the optical gap is closed, 
the exciton spectrum peak continues to exist as a reason- 
ably well-defined excitation at a slightly higher binding 
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energy. 

The exciton energy we find is not consistent with the 
0.5 eV scale proposed in Ref. |H but is reasonably con- 
sistent with the discussion in Refs. liollilllislH 

V. OPTICAL CONDUCTIVITY 

In this section we discuss the optical conductivities, 
in order to determine which features in the optical spec- 
trum may relate to the d orbitals of interest here. Also, 
previous calculation o 32 ' 33 based on the three-band model 
revealed a strong discrepancy between theory and ex- 
periment, with the theoretically calculated conductivity 
much smaller than the measured one in the region of the 
charge-transfer gap edge. Ref. [38| argued that inclusion 
of the d 3z 2_ r 2 orbital could resolve this discrepancy. 

The in-plane optical conductivities can be calculated 



from 65 

to\ - 2 JL r f rf 2 fc f(u)-f(uj + n) 

a( ' he, ]_„ n J (2tt) 2 ft 

x Tr [j(fe)A(w + ft, k)j(k)A(u, k)\ , (18) 

where Cq is the c-axis lattice parameter, /(w) is the Fermi 
function, the fe-integral is over the full Brillouin zone with 
k scaled to n divided by the in-plane lattice parameter a, 
and A(w, k) = [G(cj, k) — (w, fc)] / (2i) is the electron 
spectral function, a matrix in orbital space. We have 
chosen our Fourier transform so that the in-plane cur- 
rent operator is simply a fe-derivative of the Hamiltonian 
matrix, j = dH/dk x r^ the extra terms discussed by Ref. 
IH are not needed. 

The c-axis optical conductivity can be calculated in a 
similar manner using the current operator 



Jc = 
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and in La2Cu04 co/a ~ 1.3. 

Fig. [6] shows the calculated in-plane optical conductiv- 
ity for the T and T' phases. In the two calculations the 
A values are chosen to be at approximately the same dis- 
tance from the insulating phase boundary A C 2 so the gap 
sizes are quite similar. The two systems give very similar 
in-plane conductivities. Further, the results are very sim- 
ilar to those obtained from the three-band mode l 32 i 33 : an 
onset of absorption above around 1.8 eV and a strong ab- 
sorption at energy between 6 eV and 8 eV. The rise in the 

3 eV-5 eV range is due to the effect of the d 3z 2_ r 2 orbital. 
The introduction of additional orbitals does not increase 
the conductivity magnitude in the frequency range im- 
mediately above the gap (ft ~ 2 — 3eV) significantly: 
the disagreement with experiment previously noted in 
the three-band mode l 32 ' 33 remains. These results do not 
agree with results presented in Ref. [H. 

Fig. [7] shows the results of the c-axis optical conduc- 
tivity calculated in the T-phase. The overall magnitude 
is at least an order of magnitude smaller than that of 
the in-plane conductivity, which is a consequence of a 
much smaller hybridization strength and smaller number 
of relevant orbitals in the c-direction. In the 2 eV to 4 
eV range the conductivity is non-zero but very small. At 

4 eV the conductivity start to rise, signalling the onset 
of transitions involving the d 3z 2_ r 2 orbitals. Between 6 
eV and 8 eV there are two strong peaks which we con- 



sider to be the transition between the upper Hubbard 
band (which has a double peak structure) and the non- 
bonding p z band. 

VI. ORBITAL CHARACTER OF DOPED 
HOLES AND SHAPE OF FERMI SURFACES 

In this section we consider the variation with doping 
of the orbital character of the low-lying states. This sec- 
tion is motivated by the possibility that above a critical 
doping the d% z 2_ r 2 band begins to become occupied. 

Fig. [5] shows the doping dependence of the number 
of holes on each orbital per unit cell. From the spec- 
tral functions shown in Fig. [3] one would expect that the 
number of holes on the d 3z 2_ r 2 orbital will dramatically 
increase when the chemical potential is reduced below a 
certain point. The theoretical arguments of Ref. also 
suggest that this will occur. However the actual DMFT 
calculation is inconsistent with the rigid band picture. 
We see that while the total number of holes in the d- 
band increases linearly with doping, the hole density on 
the d 3z 2_ r 2 /p z complex remains very small even at very 
high doping levels. Thus the spectra deform as the chem- 
ical potential is reduced, in such a way that the d 3z 2_ r 2 
orbital remains far below the Fermi level. 

To gain further insight into the doping dependence we 
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FIG. 6: In-plane optical conductivities calculated for T- and 
T'-phase in the undoped case at A values with similar dis- 
tances below the insulating boundary A C 2. (7 = 9 eV and 
T = 0.1 eV. Parameters: T- phase: J = (red solid line), 
A = 21.6 eV, e d = e dz = -26.3 eV, e p = e Pz = -4.7 eV. 
J = 0.5 eV (blue dashed line), A = 19.2 eV, e d = e dz = -23.8 
eV, e p = e Pz = —4.6 eV. J — 1 eV (magenta dash-dotted 
line), A = 16.8 eV, e d = s dz = -21.3 eV, e p = e Pz = -4.5 
eV. T'-phase (black dotted line): J = 0, A = 21.8 eV, 
£d = £d z = -26.2 eV, e p = -4.4 eV. 
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FIG. 7: c-axis optical conductivities calculated for the T'- 
phase in the undoped case at A values with similar distances 
below the insulating boundary A C 2- Parameters are the same 
as in Fig. [6] 



plot the Fermi surfaces of the six-band model in T-phase 
in Fig. IH1 Panels (a) and (b) shows results obtained for 
parameters such that at half filling the model is in its 
paramagnetic insulating phase while panels (c) and (d) 
show results obtained for parameters such that at half 
filling the model is in the paramagnetic metallic phase. 



FIG. 8: Number of holes on each orbital as a function of 
the total hole doping per unit cell, calculated for the six- 
band model in T-phase. Upper panel: A = 21.6 eV (< A C 2). 
Lower panel: A = 22.8 eV (= A c i). Note that the number 
of holes on p x and p y orbital of the in-plane oxygen sites are 
combined as p x + p y (shown as blue squares with long dashed 
lines), and the number of holes on the p z orbital of the above- 
and below-plane apical oxygen sites are similarly combined 
as 2 *p z (shown as green diamonds with short dashed lines). 
Therefore at a given doping value the sum over the value at 
the four curves gives the correct total hole doping per unit 
cell. 




k,Jn kjn 



FIG. 9: Fermi Surfaces of the six-band model in T-phase. 
Panel (a): A = 21.6 eV (< A c2 ), hole doping x = 0.10. Panel 
(b): A = 21.6 eV, x = 0.35. Panel (c): A = 22.8 eV (= Ad), 
x = 0.08. Panel (d): A = 22.8 eV, x = 0.32. Parameters: 
U = 9 eV, T = 0.1 eV. Panel (a): e d = e dz = -25.3 eV, e p = 
e Pz = —3.7 eV. Panel (b): e d = e dz = —24.5 eV, e p = e Pz = 
—2.9 eV. Panel (c): e d = e dz = —26.3 eV, e p = e Pz = —3.5 
eV. Panel (d): e d = e dz = -25.5 eV, e p = e Pz = -2.7 eV. 
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hole doping 

FIG. 10: Number of holes on each orbital as a function of 
the total hole doping per unit cell, calculated for the six-band 
model in T'-phase. Upper panel: A = 21.8 eV (< A C 2). 
Lower panel: A = 22.8 eV (= A c i). Note that the number 
of holes on p x and p y orbital of the in-plane oxygen sites are 
combined as p x + p y (shown as blue squares with long dashed 
lines). 




FIG. 11: Fermi Surfaces of the six-band model in T'-phase. 
Panel (a): A = 21.8 eV (< A c2 ), hole doping x = 0.13. Panel 
(b): A = 21.8 eV, x = 0.37. Panel (c): A = 22.8 eV (= A c i), 
x = 0.12. Panel (d): A = 22.8 eV, x = 0.35. Parameters: 
U = 9 eV, T = 0.1 eV. Panel (a): e d = e dz = -25.3 eV, 
e p = -3.5 eV. Panel (b): e d = e dz = -24.5 eV, e p = -2.7 
eV. Panel (c): e d = e dz = —26.1 eV, e v = —3.3 eV. Panel 
(d): e d = e dz = -25.3 eV, e v = -2.5 eV. 



The hole doping values of panels (a) and (c) are selected 
around 0.1 and panels (b) and (d) around 0.35. We see 
that the Fermi surface includes only one sheet in all cases, 
consistent with the discussion above that the crossing of 
the Fermi energy into the d 3z 2_ r 2 band is avoided. For 
the smaller doping value the Fermi surface is approxi- 
mately a circle enclosing (tt, it) and for the larger doping 
the Fermi surface is star shaped enclosing the (0, 0) point. 
Thus, in disagreement with early slave boson studies^ 
we find that in the six-band model there is no physi- 
cally relevant doping at which holes occupy d 3z 2_ r 2/p z 
orbitals as separate bands, and the Fermi surface remains 
one-sheeted. However, we do note that the van Hove sin- 
gularity happens at around hole doping value x ~ 0.28, 
an intermediate value between what shown in panels (a) , 
(c) and (b), (d) in Fig.H 

We have repeated the entire analysis for the T" struc- 
ture, finding very similar results but with even smaller 
occupancy of the d 3z 2_ r 2 orbitals. This is understand- 
able as the hybridization to the d 3z 2_ r 2 orbital is much 
weaker once p z orbitals are removed. We also plot the 
Fermi Surface in Fig. [TT] As for the T'-phase, the Fermi 
surface has only one sheet, and the d 3z 2_ r 2 orbitals are 
not populated as separate bands. We note that the van 
Hove singularity happens at around hole doping value 
x « 0.37 which is approximately a 0.1 shift in doping 
compared to the T phase. 



VII. CONCLUSION 

In this paper, we have employed the single-site DMFT 
method to study a six-band model, which includes copper 
d x 2_ y 2 1 d 3z 2_ r 2 1 in- plane oxygen p x ,y, and (in T phase 
structure) the apical oxygen p z orbitals. This model 
is more chemically realistic than the three-band or one- 
band models frequently considered. We considered two 
structures: the T-phase, appropriate to La2CuC>4, and 
the T'-phase, appropriate to the infinite-layer cuprates 
and to the electron-doped materials such as Nd2CuC>4. 
From the model point of view these structures differ- 
ent in whether or not apical oxygen p z orbitals are in- 
duced. We have mapped out the metal/charge-transfcr- 
insulator phase diagram, finding that after the atomic- 
limit Hartree shift is accounted for, the phase bound- 
aries are systematically shifted to the metallic regime 
compared to that of the three-band model. Thus we con- 
clude that incorporating the d 3z 2_ r 2 orbital expands the 
insulating regime of the system. 

The spectral functions are calculated by analytic con- 
tinuation. The d x 2_ y 2 and p x>y spectra are observed to be 
similar to that of the three-band model. In the T-phase 
the non-bonding p z band appears as a <5-function and 
two side-bands corresponds to the bonding apical oxy- 
gen bands. Hybridization to these orbitals means that 
the d 3z 2_ r 2 has a double-peak structure. In contrast, in 
the T'-phase, the spectrum of d 3z 2_ r 2 orbital has a sin- 
gle peak. In the insulating regime, we have found that 
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the insulating gap is generically smaller in the six-band 
model than in the three-band model, for comparable cor- 
relation parameters. 

We have calculated the d-d exciton spectrum, finding 
a sharp exciton line which should be visible in experi- 
ments. In the insulating phase, the exciton states are 
inside the charge-transfer gap. In the metallic phase, 
the exciton states are at slightly higher binding energy, 
but although they overlap in energy with the tails of the 
Hubbard bands, the broadening is small. 

Both in-plane and c-axis optical conductivy are cal- 
culated. We have found, in disagreement with previ- 
ous publication, 3 - that inclusion of the additional d 3z 2_ r 2 
and apical oxygen bands does not fix the problem of the 
near-gap magnitude. The c-axis conductivity is much 
weaker and the absorption is very small in the frequency 
range 2 eV to 4 eV. Above 4 eV there is a relatively no- 
ticeable absorption due to transition from the decoupled 
apical oxygen bands to the upper Hubbard band. 

We have studied the distribution of doped holes onto 
different orbitals. We have shown that under no physi- 
cally relevant doping values that the d 3z 2_ r 2 orbitals (and 
p z orbitals in the T-phase) are populated as a separate 
band. The Fermi surface only contains one sheet. This 
is in disagreement with previous slave boson studies4I 
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Appendix 

In this appendix we present the "six-band+s" 
model involving an additional Cu 4s orbital in 



each unit cell. 



We take the basis as 



{d\\ k ,d zk ,Sk,p xk ,Pyk,Pt h k ove ,P b z t ow )- Then the Hamilto- 
nian is a 7 x 7 matrix, which may be seperated to Cu and 
O parts as 



H 



6band+s — 



H Cu 

/ H hyb y 



H 



hyb 

6band+s 



H° 

AJ -6band+s 



(A-l) 



where 



I s d 

H6band+s = " £ d z 

\ e 

f£ P 

e p 

e Pz 

V 



ttO 

6band+s 




(A-2) 
(A-3) 



and the hybridization between Cu and O orbitals: 



H 



hyb 

6band+s 

2it p d sin 4* — 2it p d sin -y- 

-2it pdz sin -2it pdz sin -f 



2it ps sin 4f 



2it ps sin -f 





t PzS 







(A-4) 



The downfolding 50 of Eq. (|A-1|) leads to the Hamilto- 
nian considered in the main text. The effective oxygen- 
oxygen hopping amplitudes are 



ps 



^pp 



b PzPz 



e s — ef 

tpS ' tp z S 



Pz S 



e F 



Note that this implies that t pPz /t pp 



tp z p z I^PPz 



(A-5) 
(A-6) 

(A-7) 
which 



has been used in obtaining the value of t PzPz in the main 
text. 
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